ABSTRACT
Background: The time-varying reproduction number (Rt) is an important measure of epidemic transmissibility; it can directly inform policy decisions and the optimisation of control measures. EpiEstim is a widely used software tool that uses case incidence and the serial interval (SI, time between symptoms in a case and their infector) to estimate Rt in real-time. The incidence and the SI distribution must be provided at the same temporal resolution, which limits the applicability of EpiEstim and other similar methods, e.g. for pathogens with a mean SI shorter than the frequency of incidence reporting. Methods We use an expectation-maximisation algorithm to reconstruct daily incidence from temporally aggregated data, from which Rt can then be estimated using EpiEstim. We assess the validity of our method using an extensive simulation study and apply it to COVID-19 and influenza data. The method is implemented in the opensource R package EpiEstim. Findings For all datasets, the influence of intra-weekly variability in reported data was mitigated by using aggregated weekly data. Rt estimated on weekly sliding windows using incidence reconstructed from weekly data was strongly correlated with estimates from the original daily data. The simulation study revealed that Rt was well estimated in all scenarios and regardless of the temporal aggregation of the data. In the presence of weekend effects, Rt estimates from reconstructed data were more successful at recovering the true value of Rt than those obtained from reported daily data. Interpretation Rt can be successfully recovered from aggregated data, and estimation accuracy can even be improved by smoothing out administrative noise in the reported data. Funding MRC doctoral training partnership, MRC centre for global infectious disease analysis, the NIHR HPRU in Modelling and Health Economics, and the Academy of Medical Sciences Springboard, funded by the AMS, Wellcome Trust, BEIS, the British Heart Foundation and Diabetes UK.
Subject(s)
COVID-19 , Diabetes MellitusABSTRACT
Recent months have demonstrated that emerging variants may set back the global COVID-19 response. The ability to rapidly assess the threat of new variants in real-time is critical for timely optimisation of control strategies. We extend the EpiEstim R package, designed to estimate the time-varying reproduction number (Rt), to estimate in real-time the effective transmission advantage of a new variant compared to a reference variant. Our method can combine information across multiple locations and over time and was validated using an extensive simulation study, designed to mimic a variety of real-time epidemic contexts. We estimate that the SARS-CoV-2 Alpha variant is 1.46 (95% Credible Interval 1.44-1.47) and 1.29, (95% CrI 1.29-1.30) times more transmissible than the wild type, using data from England and France respectively. We further estimate that Beta and Gamma combined are 1.25 (95% CrI 1.24-1.27) times more transmissible than the wildtype (France data). All results are in line with previous estimates from literature, but could have been obtained earlier and more easily with our off-the-shelf open-source tool. Our tool can be used as an important first step towards quantifying the threat of new variants in real-time. Given the popularity of EpiEstim, this extension will likely be used widely to monitor the co-circulation and/or emergence of multiple variants of infectious pathogens.
Subject(s)
COVID-19ABSTRACT
We report on the second and final part of a pre-registered forecasting study on COVID-19 cases and deaths in Germany and Poland. Fifteen independent research teams provided forecasts at lead times of one through four weeks from January through mid-April 2021. Compared to the first part (October--December 2020), the number of participating teams increased, and a number of teams started providing subnational-level forecasts. The addressed time period is characterized by rather stable non-pharmaceutical interventions in both countries, making short-term predictions more straightforward than in the first part of our study. In both countries, case counts declined initially, before rebounding due to the rise of the B.1.1.7 variant. Deaths declined through most of the study period in Germany while in Poland they increased after a prolonged plateau. Many, though not all, models outperformed a simple baseline model up to four weeks ahead, with ensemble methods showing very good relative performance. Major trend changes in reported cases, however, remained challenging to predict.
Subject(s)
COVID-19ABSTRACT
From 8th March to 29th November 2020, we produced weekly estimates of SARS-CoV-2 transmissibility and forecasts of deaths due to COVID-19 for 81 countries with evidence of sustained transmission. We also developed a novel heuristic to combine weekly estimates of transmissibility to produce forecasts over a 4-week horizon. We evaluated the robustness of the forecasts using relative error, coverage probability, and comparisons with null models. During the 39-week period covered by this study, both the short- and medium-term forecasts captured well the epidemic trajectory across different waves of COVID-19 infections with small relative errors over the forecast horizon. The model was well calibrated with 56.3\% and 45.6\% of the observations lying in the 50\% Credible Interval in 1-week and 4-week ahead forecasts respectively. We could accurately characterise the overall phase of the epidemic up to 4-weeks ahead in 84.9\% of country-days. The medium-term forecasts can be used in conjunction with the short-term forecasts of COVID-19 mortality as a useful planning tool as countries continue to relax public health measures.
Subject(s)
COVID-19ABSTRACT
BackgroundAs of July 2021, more than 180,000,000 cases of COVID-19 have been reported across the world, with more than 4 million deaths. Mathematical modelling and forecasting efforts have been widely used to inform policy-making and to create situational awareness. Methods and FindingsFrom 8th March to 29th November 2020, we produced weekly estimates of SARS-CoV-2 transmissibility and forecasts of deaths due to COVID-19 for countries with evidence of sustained transmission. The estimates and forecasts were based on an ensemble model comprising of three models that were calibrated using only the reported number of COVID-19 cases and deaths in each country. We also developed a novel heuristic to combine weekly estimates of transmissibility and potential changes in population immunity due to infection to produce forecasts over a 4-week horizon. We evaluated the robustness of the forecasts using relative error, coverage probability, and comparisons with null models. ConclusionsDuring the 39-week period covered by this study, we produced short- and medium-term forecasts for 81 countries. Both the short- and medium-term forecasts captured well the epidemic trajectory across different waves of COVID-19 infections with small relative errors over the forecast horizon. The model was well calibrated with 56.3% and 45.6% of the observations lying in the 50% Credible Interval in 1-week and 4-week ahead forecasts respectively. We could accurately characterise the overall phase of the epidemic up to 4-weeks ahead in 84.9% of country-days. The medium-term forecasts can be used in conjunction with the short-term forecasts of COVID-19 mortality as a useful planning tool as countries continue to relax stringent public health measures that were implemented to contain the pandemic.
Subject(s)
COVID-19 , DeathABSTRACT
Measuring COVID-19 spread remains challenging in many countries due to testing limitations. In Java, reported cases and deaths increased throughout 2020 despite intensive control measures, particularly within Jakarta and during Ramadan. However, underlying trends are likely obscured by variations in case ascertainment. COVID-19 protocol funerals in Jakarta provide alternative data indicating a substantially higher burden than observed within confirmed deaths. Transmission estimates using this metric follow mobility trends, suggesting earlier and more sustained intervention impact than observed in routine data. Modelling suggests interventions have lessened spread to rural, older communities with weaker healthcare systems, though predict healthcare capacity will soon be exceeded in much of Java without further control. Our results highlight the important role syndrome-based measures of mortality can play in understanding COVID-19 transmission and burden.
Subject(s)
COVID-19ABSTRACT
As of 1st June 2020, the US Centers for Disease Control and Prevention reported 104,232 confirmed or probable COVID-19-related deaths in the US. This was more than twice the number of deaths reported in the next most severely impacted country. We jointly modelled the US epidemic at the state-level, using publicly available death data within a Bayesian hierarchical semi-mechanistic framework. For each state, we estimate the number of individuals that have been infected, the number of individuals that are currently infectious and the time-varying reproduction number (the average number of secondary infections caused by an infected person). We used changes in mobility to capture the impact that non-pharmaceutical interventions and other behaviour changes have on the rate of transmission of SARS-CoV-2. Nationally, we estimated 3.7% [3.4%-4.0%] of the population had been infected by 1st June 2020, with wide variation between states, and approximately 0.01% of the population was infectious. We also demonstrated that good model forecasts of deaths for the next 3 weeks with low error and good coverage of our credible intervals.
Subject(s)
COVID-19 , Coinfection , Oculocerebrorenal Syndrome , DeathABSTRACT
Brazil is an epicentre for COVID-19 in Latin America. In this report we describe the Brazilian epidemic using three epidemiological measures: the number of infections, the number of deaths and the reproduction number. Our modelling framework requires sufficient death data to estimate trends, and we therefore limit our analysis to 16 states that have experienced a total of more than fifty deaths. The distribution of deaths among states is highly heterogeneous, with 5 states---Sao Paulo, Rio de Janeiro, Ceara, Pernambuco and Amazonas---accounting for 81% of deaths reported to date. In these states, we estimate that the percentage of people that have been infected with SARS-CoV-2 ranges from 3.3% (95% CI: 2.8%-3.7%) in Sao Paulo to 10.6% (95% CI: 8.8%-12.1%) in Amazonas. The reproduction number (a measure of transmission intensity) at the start of the epidemic meant that an infected individual would infect three or four others on average. Following non-pharmaceutical interventions such as school closures and decreases in population mobility, we show that the reproduction number has dropped substantially in each state. However, for all 16 states we study, we estimate with high confidence that the reproduction number remains above 1. A reproduction number above 1 means that the epidemic is not yet controlled and will continue to grow. These trends are in stark contrast to other major COVID-19 epidemics in Europe and Asia where enforced lockdowns have successfully driven the reproduction number below 1. While the Brazilian epidemic is still relatively nascent on a national scale, our results suggest that further action is needed to limit spread and prevent health system overload.
Subject(s)
COVID-19 , Death , InfectionsABSTRACT
Italy was the first European country to experience sustained local transmission of COVID-19. As of 1st May 2020, the Italian health authorities reported 28,238 deaths nationally. To control the epidemic, the Italian government implemented a suite of non-pharmaceutical interventions (NPIs), including school and university closures, social distancing and full lockdown involving banning of public gatherings and non essential movement. In this report, we model the effect of NPIs on transmission using data on average mobility. We estimate that the average reproduction number (a measure of transmission intensity) is currently below one for all Italian regions, and significantly so for the majority of the regions. Despite the large number of deaths, the proportion of population that has been infected by SARS-CoV-2 (the attack rate) is far from the herd immunity threshold in all Italian regions, with the highest attack rate observed in Lombardy (13.18% [10.66%-16.70%]). Italy is set to relax the currently implemented NPIs from 4th May 2020. Given the control achieved by NPIs, we consider three scenarios for the next 8 weeks: a scenario in which mobility remains the same as during the lockdown, a scenario in which mobility returns to pre-lockdown levels by 20%, and a scenario in which mobility returns to pre-lockdown levels by 40%. The scenarios explored assume that mobility is scaled evenly across all dimensions, that behaviour stays the same as before NPIs were implemented, that no pharmaceutical interventions are introduced, and it does not include transmission reduction from contact tracing, testing and the isolation of confirmed or suspected cases. New interventions, such as enhanced testing and contact tracing are going to be introduced and will likely contribute to reductions in transmission; therefore our estimates should be viewed as pessimistic projections. We find that, in the absence of additional interventions, even a 20% return to pre-lockdown mobility could lead to a resurgence in the number of deaths far greater than experienced in the current wave in several regions. Future increases in the number of deaths will lag behind the increase in transmission intensity and so a second wave will not be immediately apparent from just monitoring of the daily number of deaths. Our results suggest that SARS-CoV-2 transmission as well as mobility should be closely monitored in the next weeks and months. To compensate for the increase in mobility that will occur due to the relaxation of the currently implemented NPIs, adherence to the recommended social distancing measures alongside enhanced community surveillance including swab testing, contact tracing and the early isolation of infections are of paramount importance to reduce the risk of resurgence in transmission.